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I 1 Abstract 

Recent theoretical and empirical work in statistical machine learning has demonstrated 
, ^ , the importance of learning algorithms for deep architectures, i.e., function classes obtained 

by composing multiple non-linear transformations. Self-taught learning (exploiting unla- 
i-H beled examples or examples from other distributions) has already been applied to deep 

learners, but mostly to show the advantage of unlabeled examples. Here we explore the 
ON advantage brought by out-of-distribution examples. For this purpose we developed a pow- 

OO erful generator of stochastic variations and noise processes for character images, including 

not only affine transformations but also slant, local elastic deformations, changes in thick- 
ness, background images, grey level changes, contrast, occlusion, and various types of 
ON noise. The out-of-distribution examples are obtained from these highly distorted images 

or by including examples of object classes different from those in the target test set. We 
show that deep learners benefit more from out-of-distribution examples than a correspond- 
ing shallow learner, at least in the area of handwritten character recognition. In fact, we 
show that they beat previously published results and reach human-level performance on 
both handwritten digit classification and 62-class handwritten character recognition. 
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1 Introduction 



Deep Learning has emerged as a promising new area of research in statistical machine learn- 
ing (see Bengio [1 J for a review). Learning algorithms for deep architectures are centered on 
the learning of useful representations of data, which are better suited to the task at hand, and 
are organized in a hierarchy with multiple levels. This is in part inspired by observations of 
the mammalian visual cortex, which consists of a chain of processing elements, each of which 
is associated with a different representation of the raw visual input. In fact, it was found re- 
cently that the features learnt in deep architectures resemble those observed in the first two of 
these stages (in areas VI and V2 of visual cortex) Q, and that they become more and more 
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invariant to factors of variation (such as camera movement) in higher layers [3]. Learning a 
hierarchy of features increases the ease and practicality of developing representations that are 
at once tailored to specific tasks, yet are able to borrow statistical strength from other related 
tasks (e.g., modeling different kinds of objects). Finally, learning the feature representation 
can lead to higher-level (more abstract, more general) features that are more robust to unan- 
ticipated sources of variance extant in real data. 

Self-taught learning |4| is a paradigm that combines principles of semi-supervised and 
multi-task learning: the learner can exploit examples that are unlabeled and possibly come 
from a distribution different from the target distribution, e.g., from other classes than those of 
interest. It has already been shown that deep learners can clearly take advantage of unsuper- 
vised learning and unlabeled examples [ 1 1 5 1, but more needs to be done to explore the impact 
of out-of -distribution examples and of the multi-task setting (one exception is |6|, which uses 
a different kind of learning algorithm). In particular the relative advantage of deep learning 
for these settings has not been evaluated. The hypothesis discussed in the conclusion is that a 
deep hierarchy of features may be better able to provide sharing of statistical strength between 
different regions in input space or different tasks. 

Whereas a deep architecture can in principle be more powerful than a shallow one in terms 
of representation, depth appears to render the training problem more difficult in terms of opti- 
mization and local minima. It is also only recently that successful algorithms were proposed 
to overcome some of these difficulties. All are based on unsupervised learning, often in an 
greedy layer- wise "unsupervised pre-training" stage [ 1 ]. One of these layer initialization tech- 
niques, applied here, is the Denoising Auto-encoder (DA) Q (see Figure[T|, which performed 
similarly or better than previously proposed Restricted Boltzmann Machines in terms of un- 
supervised extraction of a hierarchy of features useful for classification. Each layer is trained 
to denoise its input, creating a layer of features that can be used as input for the next layer. 

In this paper we ask the following questions: 

• Do the good results previously obtained with deep architectures on the MNIST digit 
images generalize to the setting of a much larger and richer (but similar) dataset, the NIST 
special database 19, with 62 classes and around 800k examples? 

• To what extent does the perturbation of input images (e.g. adding noise, affine trans- 
formations, background images) make the resulting classifiers better not only on similarly 
perturbed images but also on the original clean examples'] We study this question in the 
context of the 62-class and 10-class tasks of the NIST special database 19. 

• Do deep architectures benefit more from such out-of-distribution examples, i.e. do they 
benefit more from the self-taught learning JU framework? We use highly perturbed examples 
to generate out-of-distribution examples. 

• Similarly, does the feature learning step in deep learning algorithms benefit more from 
training with moderately different classes (i.e. a multi-task learning scenario) than a corre- 
sponding shallow and purely supervised architecture? We train on 62 classes and test on 10 
(digits) or 26 (upper case or lower case) to answer this question. 

Our experimental results provide positive evidence towards all of these questions, as well 
as classifiers that reach human-level performance on 62-class isolated character recognition 
and beat previously published results on the NIST dataset (special database 19). To achieve 
these results, we introduce in the next section a sophisticated system for stochastically trans- 
forming character images and then explain the methodology, which is based on training with 
or without these transformed images and testing on clean ones. We measure the relative ad- 
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vantage of out-of-distribution examples (perturbed or out-of-class) for a deep learner vs a 
supervised shallow one. Code for generating these transformations as well as for the deep 
learning algorithms are made available at |http : //hg . assembla . com/ if 1 62 6~6] We 
estimate the relative advantage for deep learners of training with other classes than those of 
interest, by comparing learners trained with 62 classes with learners trained with only a subset 
(on which they are then tested). The conclusion discusses the more general question of why 
deep learners may benefit so much from the self-taught learning framework. Since out-of- 
distribution data (perturbed or from other related classes) is very common, this conclusion is 
of practical importance. 



2 Perturbed and Transformed Character Images 
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This section describes the different transformations we used to stochastically 
transform 32 x 32 source images (such as the one on the left) in order to 
obtain data from a larger distribution which covers a domain substantially 
larger than the clean characters distribution from which we start. Although 
character transformations have been used before to improve character rec- 
ognizers, this effort is on a large scale both in number of classes and in 
the complexity of the transformations, hence in the complexity of the leam- 
Original tas j f CQ( j e f or ^ggg transformations (mostly python) is available at 

|http : / /hg . assembla . com/ if 1 62 6 6| All the modules in the pipeline share a global 
control parameter (0 < complexity < 1) that allows one to modulate the amount of deforma- 
tion or noise introduced. There are two main parts in the pipeline. The first one, from slant to 
pinch below, performs transformations. The second part, from blur to contrast, adds different 
kinds of noise. 



2.1 Transformations 



Thickness 



To change character thickness, morphological operators of dilation and ero- 
sion |8] [9) are applied. The neighborhood of each pixel is multiplied element- 
wise with a structuring element matrix. The pixel value is replaced by the max- 
imum or the minimum of the resulting matrix, respectively for dilation or ero- 
sion. Ten different structural elements with increasing dimensions (largest is 
5x5) were used. For each image, randomly sample the operator type (dilation 
or erosion) with equal probability and one structural element from a subset of the 
n = round(m x complexity) smallest structuring elements where m = 10 for 
dilation and m = 6 for erosion (to avoid completely erasing thin characters). A 
neutral element (no transformation) is always present in the set. 



3 



Slant 



To produce slant, each row of the image is shifted proportionally to its height: 

^ shift = round(slant x height), slant ~ U[— complexity , complexity]. 
The shift is randomly chosen to be either to the left or to the right. 

Affine Transformations 

A 2 x 3 affine transform matrix (with parameters (a, b, c, d, e, /)) is sampled 
according to the complexity. Output pixel (x, y) takes the value of input 

E^H pixel nearest to (ax + by + c, dx + ey + /), producing scaling, transla- 
tion, rotation and shearing. Marginal distributions of (a,b,c,d,e, /) have 
been tuned to forbid large rotations (to avoid confusing classes) but to give 
good variability of the transformation: a and d ~ t/[l — ^complexity , 1 + 
3 complexity], b and e ~ {/[— 3 complexity, 3 complexity], and c and 
/ ~ [/[— 4 complexity, 4 complexity]. 



Local Elastic Deformations 



The local elastic deformation module induces a "wiggly" effect in the image, 
following Simard et al. ifTUl . which provides more details. The intensity of the 
displacement fields is given by a = {/complexity x 10.0, which are convolved 
with a Gaussian 2D kernel (resulting in a blur) of standard deviation a = 10 — 

7 x {/complexity. 



Pinch 
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The pinch module applies the "Whirl and pinch" GIMP filter with whirl set to 0. 
A pinch is "similar to projecting the image onto an elastic surface and pressing or 
pulling on the center of the surface" (GIMP documentation manual). For a square 
input image, draw a radius-r disk around its center C. Any pixel P belonging 
to that disk has its value replaced by the value of a "source" pixel in the original 
image, on the line that goes through C and P, but at some other distance cfe. 
Define di = distance(P,C) and d 2 = sin( ! ^ x )~ pinch x d±, where pinch 
is a parameter of the filter. The actual value is given by bilinear interpolation 
considering the pixels around the (non-integer) source position thus found. Here 
pinch ~ U[— complexity, 0.7 x complexity]. 
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2.2 Injecting Noise 
Motion Blur 



i 
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Occlusion 
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The motion blur module is GIMP's "linear motion blur", which has parame- 
ters length and angle. The value of a pixel in the final image is approximately 
the mean of the first length pixels found by moving in the angle direction, 
angle ~ U[0, 360] degrees, and length ~ Normal(0, (3 x complexity) 2 ). 



The occlusion module selects a random rectangle from an occluder character 
image and places it over the original occluded image. Pixels are combined 
by taking the max(occluder, occluded), i.e. keeping the lighter ones. The 
rectangle corners are sampled so that larger complexity gives larger rectangles. 
The destination position in the occluded image are also sampled according to 
a normal distribution. This module is skipped with probability 60%. 



Gaussian Smoothing 

With the Gaussian smoothing module, different regions of the image are 
spatially smoothed. This is achieved by first convolving the image with an 
isotropic Gaussian kernel of size and variance chosen uniformly in the ranges 

H[12, 12 + 20 x complexity] and [2, 2 + 6 x complexity]. This filtered image 
is normalized between and 1. We also create an isotropic weighted averaging 
window, of the kernel size, with maximum value at the center. For each image 
we sample uniformly from 3 to 3 + 10 x complexity pixels that will be averaging 
centers between the original image and the filtered one. We initialize to zero a 
mask matrix of the image size. For each selected pixel we add to the mask the 
averaging window centered on it. The final image is computed from the follow- 
ing element-wise operation: ' ma9e+/ '"^ a9exmgst ■ This module is skipped 
with probability 75%. 

Permute Pixels 



This module permutes neighbouring pixels. It first selects a fraction com P l ^ xlt v 
of pixels randomly in the image. Each of these pixels is then sequentially 
exchanged with a random pixel among its four nearest neighbors (on its left, 
right, top or bottom). This module is skipped with probability 80%. 
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Gaussian Noise 
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Background 


I 
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The Gaussian noise module simply adds, to each pixel of the image indepen- 
dently, a noise - Normal{0, (£2^£^*a)2) This moc jule is skipped with prob- 
ability 70%. 



Image Addition 

Following Larochelle et al. ifTTI . the background image module adds a ran- 
dom background image behind the letter, from a randomly chosen natural im- 
age, with contrast adjustments depending on complexity, to preserve more or 
less of the original character image. 



Salt and Pepper Noise 

The salt and pepper noise module adds noise ~ U[0, 1] to random subsets of 
pixels. The number of selected pixels is 0.2 x complexity. This module is 
skipped with probability 75%. 



Scratches 

The scratches module places line-like white patches on the image. The lines are 
heavily transformed images of the digit "1" (one), chosen at random among 500 
such 1 images, randomly cropped and rotated by an angle ~ Normal(0, (100 x 
complexity) 2 (in degrees), using bi-cubic interpolation. Two passes of a grey- 
scale morphological erosion filter are applied, reducing the width of the line by 
an amount controlled by complexity. This module is skipped with probability 
85%. The probabilities of applying 1, 2, or 3 patches are (50%,30%,20%). 

Grey Level and Contrast Changes 





The grey level and contrast module changes the contrast by changing grey 
levels, and may invert the image polarity (white to black and black to white). 
The contrast is C ~ U[l — 0.85 x complexity, 1] so the image is normalized 
into [^ttS 1 — ^-ip-}- The polarity is inverted with probability 50%. 



3 Experimental Setup 

Much previous work on deep learning had been performed on the MNIST digits task lfT2l[T3l 
[T4l[T5ll . with 60 000 examples, and variants involving 10 000 examples ifTBIfTTIl . The focus 
here is on much larger training sets, from 10 times to to 1000 times larger, and 62 classes. 
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The first step in constructing the larger datasets (called NISTP and P07) is to sample from 
a data source: NIST (NIST database 19), Fonts, Captchas, and OCR data (scanned machine 
printed characters). Once a character is sampled from one of these sources (chosen randomly), 
the second step is to apply a pipeline of transformations and/or noise processes described in 
section |2 

To provide a baseline of error rate comparison we also estimate human performance on 
both the 62-class task and the 10-class digits task. We compare the best Multi-Layer Percep- 
trons (MLP) against the best Stacked Denoising Auto-encoders (SDA), when both models' 
hyper-parameters are selected to minimize the validation set error. We also provide a compar- 
ison against a precise estimate of human performance obtained via Amazon's Mechanical Turk 
(AMT) service (http://mturk.com). AMT users are paid small amounts of money to perform 
tasks for which human intelligence is required. Mechanical Turk has been used extensively 
in natural language processing and vision. AMT users were presented with 10 character im- 
ages (from a test set) and asked to choose 10 corresponding ASCII characters. They were 
forced to choose a single character class (either among the 62 or 10 character classes) for each 
image. 80 subjects classified 2500 images per (dataset,task) pair. Different humans labelers 
sometimes provided a different label for the same example, and we were able to estimate the 
error variance due to this effect because each image was classified by 3 different persons. The 
average error of humans on the 62-class task NIST test set is 18.2%, with a standard error of 
0.1%. 

3.1 Data Sources 

NIST. Our main source of characters is the NIST Special Database 19 lfl8l . widely used for 
training and testing character recognition systems [19, 201I2T1I221 . The dataset is composed 
of 8 14255 digits and characters (upper and lower cases), with hand checked classifications, 
extracted from handwritten sample forms of 3600 writers. The characters are labelled by 
one of the 62 classes corresponding to "0"-"9" 7 "A"-"Z" and "a"-"z". The dataset contains 8 
parts (partitions) of varying complexity. The fourth partition (called hsf^, 82587 examples), 
experimentally recognized to be the most difficult one, is the one recommended by NIST 
as a testing set and is used in our work as well as some previous work lfT9l l20l I2T1 l22l for 
that purpose. We randomly split the remainder (731668 examples) into a training set and a 
validation set for model selection. The performances reported by previous work on that dataset 
mostly use only the digits. Here we use all the classes both in the training and testing phase. 
This is especially useful to estimate the effect of a multi-task setting. The distribution of the 
classes in the NIST training and test sets differs substantially, with relatively many more digits 
in the test set, and a more uniform distribution of letters in the test set (whereas in the training 
set they are distributed more like in natural text). 

Fonts. In order to have a good variety of sources we downloaded an important number 
of free fonts from: |http : / /eg . scs . car let on . ca / ~| luc/ f reef onts . html. In- 
cluding the operating system's (Windows 7) fonts, there is a total of 9817 different fonts that 
we can choose uniformly from. The chosen ttf file is either used as input of the Captcha 
generator (see next item) or, by producing a corresponding image, directly as input to our 
models. 

Captchas. The Captcha data source is an adaptation of the pycaptcha library (a python 
based captcha generator library) for generating characters of the same format as the NIST 
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dataset. This software is based on a random character class generator and various kinds of 
transformations similar to those described in the previous sections. In order to increase the 
variability of the data generated, many different fonts are used for generating the characters. 
Transformations (slant, distortions, rotation, translation) are applied to each randomly gener- 
ated character with a complexity depending on the value of the complexity parameter provided 
by the user of the data source. 

OCR data. A large set (2 million) of scanned, OCRed and manually verified machine- 
printed characters where included as an additional source. This set is part of a larger cor- 
pus being collected by the Image Understanding Pattern Recognition Research group led by 
Thomas Breuel at University of Kaiserslautern ( |http : / /www . iupr . com| , and which will 
be publicly released. 

3.2 Data Sets 

All data sets contain 32x32 grey -level images (values in [0, 1]) associated with a label from 
one of the 62 character classes. 

NIST. This is the raw NIST special database 19 0D- It has {651668 / 80000 / 82587} 
{training / validation / test} examples. 

P07. This dataset is obtained by taking raw characters from all four of the above sources 
and sending them through the transformation pipeline described in section [2] For each new 
example to generate, a data source is selected with probability 10% from the fonts, 25% from 
the captchas, 25% from the OCR data and 40% from NIST. We apply all the transformations 
in the order given above, and for each of them we sample uniformly a complexity in the range 
[0, 0.7]. It has {81920000 / 80000 / 20000} {training / validation / test} examples. 

NISTP. This one is equivalent to P07 (complexity parameter of 0.7 with the same propor- 
tions of data sources) except that we only apply transformations from slant to pinch. There- 
fore, the character is transformed but no additional noise is added to the image, giving images 
closer to the NIST dataset. It has {81920000 / 80000 / 20000} {training / validation / test} 
examples. 

3.3 Models and their Hyperparameters 

The experiments are performed using MLPs (with a single hidden layer) and SDAs. Hyper- 
parameters are selected based on the NISTP validation set error. 

Multi-Layer Perceptrons (MLP). Whereas previous work had compared deep architec- 
tures to both shallow MLPs and SVMs, we only compared to MLPs here because of the very 
large datasets used (making the use of SVMs computationally challenging because of their 
quadratic scaling behavior). Preliminary experiments on training SVMs (libSVM) with sub- 
sets of the training set allowing the program to fit in memory yielded substantially worse 
results than those obtained with MLPs. For training on nearly a billion examples (with the 
perturbed data), the MLPs and SDA are much more convenient than classifiers based on ker- 
nel methods. The MLP has a single hidden layer with tanli activation functions, and softmax 
(normalized exponentials) on the output layer for estimating P(class\image). The number of 
hidden units is taken in {300, 500, 800, 1000, 1500}. Training examples are presented in mini- 
batches of size 20. A constant learning rate was chosen among {0.001, 0.01, 0.025, 0.075, 0.1, 0.5}. 
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Stacked Denoising Auto-Encoders (SDA). Various auto-encoder variants and Restricted 
Boltzmann Machines (RBMs) can be used to initialize the weights of each layer of a deep MLP 
(with many hidden layers) lfT2l[T3lfT4l . apparently setting parameters in the basin of attraction 
of supervised gradient descent yielding better generalization ||23l . This initial unsupervised 
pre-training phase uses all of the training images but not the training labels. Each layer is 
trained in turn to produce a new representation of its input (starting from the raw pixels). It 
is hypothesized that the advantage brought by this procedure stems from a better prior, on the 
one hand taking advantage of the link between the input distribution P(x) and the conditional 
distribution of interest P(y\x) (like in semi-supervised learning), and on the other hand taking 
advantage of the expressive power and bias implicit in the deep architecture (whereby complex 
concepts are expressed as compositions of simpler ones through a deep hierarchy). 
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Figure 1 : Illustration of the computations and training criterion for the denoising auto-encoder 
used to pre-train each layer of the deep architecture. Input x of the layer (i.e. raw input or 
output of previous layer) s corrupted into x and encoded into code y by the encoder fe(-). 
The decoder ge>(-) maps y to reconstruction z, which is compared to the uncorrupted input x 
through the loss function Lh{x, z), whose expected value is approximately minimized during 
training by tuning and 6' . 

Here we chose to use the Denoising Auto-encoder [ 17 1 as the building block for these deep 
hierarchies of features, as it is simple to train and explain (see Figure [T] as well as tutorial and 
code there: [http : / /deeplearning . net /tutorial] !, provides efficient inference, and 
yielded results comparable or better than RBMs in series of experiments ifTTll . During training, 
a Denoising Auto-encoder is presented with a stochastically corrupted version of the input 
and trained to reconstruct the uncorrupted input, forcing the hidden units to represent the 
leading regularities in the data. Here we use the random binary masking corruption (which 
sets to a random subset of the inputs). Once it is trained, in a purely unsupervised way, 
its hidden units' activations can be used as inputs for training a second one, etc. After this 
unsupervised pre-training stage, the parameters are used to initialize a deep MLR which is 
fine-tuned by the same standard procedure used to train them (see previous section). The SDA 
hyper-parameters are the same as for the MLR with the addition of the amount of corruption 
noise (we used the masking noise process, whereby a fixed proportion of the input values, 
randomly selected, are zeroed), and a separate learning rate for the unsupervised pre-training 
stage (selected from the same above set). The fraction of inputs corrupted was selected among 
{10%, 20%, 50%}. Another hyper-parameter is the number of hidden layers but it was fixed 
to 3 based on previous work with SDAs on MNIST 1 17]. The size of the hidden layers was 
kept constant across hidden layers, and the best results were obtained with the largest values 
that we could experiment with given our patience, with 1000 hidden units. 
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Figure 2: SDAx are the deep models. Error bars indicate a 95% confidence interval. indi- 
cates that the model was trained on NIST, 1 on NISTP, and 2 on P07. Left: overall results of 
all models, on NIST and NISTP test sets. Right: error rates on NIST test digits only, along 
with the previous results from literature lfl9l l20l I2T1 l22l respectively based on ART, nearest 
neighbors, MLPs, and SVMs. 

4 Experimental Results 

The models are either trained on NIST (MLPO and SDAO), NISTP (MLP1 and SDA1), or P07 
(MLP2 and SDA2), and tested on either NIST, NISTP or P07, either on the 62-class task or 
on the 10-digits task. Training (including about half for unsupervised pre-training, for DAs) 
on the larger datasets takes around one day on a GPU-285. Figure [2] summarizes the results 
obtained, comparing humans, the three MLPs (MLPO, MLP1, MLP2) and the three SDAs 
(SDAO, SDA1, SDA2), along with the previous results on the digits NIST special database 
19 test set from the literature, respectively based on ARTMAP neural networks [19|, fast 
nearest-neighbor search 1 20 1, MLPs lETTl . and SVMs [22 J. More detailed and complete nu- 
merical results (figures and tables, including standard errors on the error rates) can be found 
in Appendix. The deep learner not only outperformed the shallow ones and previously pub- 
lished performance (in a statistically and qualitatively significant way) but when trained with 
perturbed data reaches human performance on both the 62-class task and the 10-class (digits) 
task. 17% error (SDA1) or 18% error (humans) may seem large but a large majority of the 
errors from humans and from SDA1 are from out-of-context confusions (e.g. a vertical bar 
can be a "1", an "1" or an "L", and a "c" and a "C" are often indistinguishible). 

In addition, as shown in the left of Figure [3] the relative improvement in error rate brought 
by self-taught learning is greater for the SDA, and these differences with the MLP are sta- 
tistically and qualitatively significant. The left side of the figure shows the improvement to 
the clean NIST test set error brought by the use of out-of-distribution examples (i.e. the per- 
turbed examples examples from NISTP or P07). Relative percent change is measured by 
taking 100% x (original model's error / perturbed-data model's error - 1). The right side of 
Figure [3] shows the relative improvement brought by the use of a multi-task setting, in which 
the same model is trained for more classes than the target classes of interest (i.e. training with 
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Figure 3: Relative improvement in error rate due to self-taught learning. Left: Improvement 
(or loss, when negative) induced by out-of-distribution examples (perturbed data). Right: 
Improvement (or loss, when negative) induced by multi-task learning (training on all classes 
and testing only on either digits, upper case, or lower-case). The deep learner (SDA) benefits 
more from both self-taught learning scenarios, compared to the shallow MLR 



all 62 classes when the target classes are respectively the digits, lower-case, or upper-case 
characters). Again, whereas the gain from the multi-task setting is marginal or negative for 
the MLR it is substantial for the SDA. Note that to simplify these multi-task experiments, only 
the original NIST dataset is used. For example, the MLP-digits bar shows the relative percent 
improvement in MLP error rate on the NIST digits test set is 100% x (single-task model's 
error / multi-task model's error - 1). The single-task model is trained with only 10 outputs 
(one per digit), seeing only digit examples, whereas the multi-task model is trained with 62 
outputs, with all 62 character classes as examples. Hence the hidden units are shared across 
all tasks. For the multi-task model, the digit error rate is measured by comparing the correct 
digit class with the output class associated with the maximum conditional probability among 
only the digit classes outputs. The setting is similar for the other two target classes (lower case 
characters and upper case characters). 



5 Conclusions and Discussion 

We have found that the self-taught learning framework is more beneficial to a deep learner 
than to a traditional shallow and purely supervised learner. More precisely, the answers are 
positive for all the questions asked in the introduction. 

• Do the good results previously obtained with deep architectures on the MNIST 
digits generalize to a much larger and richer (but similar) dataset, the NIST special 
database 19, with 62 classes and around 800k examples? Yes, the SDA systematically 
outperformed the MLP and all the previously published results on this dataset (the ones that 
we are aware of), in fact reaching human-level performance at around 17% error on the 62- 
class task and 1.4% on the digits, and beating previously published results on the same data. 

• To what extent do self-taught learning scenarios help deep learners, and do they 
help them more than shallow supervised ones? We found that distorted training examples 
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not only made the resulting classifier better on similarly perturbed images but also on the 
original clean examples, and more importantly and more novel, that deep architectures benefit 
more from such out- of- distribution examples. MLPs were helped by perturbed training exam- 
ples when tested on perturbed input images (65% relative improvement on NISTP) but only 
marginally helped (5% relative improvement on all classes) or even hurt (10% relative loss on 
digits) with respect to clean examples . On the other hand, the deep SDAs were significantly 
boosted by these out-of-distribution examples. Similarly, whereas the improvement due to 
the multi-task setting was marginal or negative for the MLP (from +5.6% to -3.6% relative 
change), it was quite significant for the SDA (from +13% to +27% relative change), which 
may be explained by the arguments below. 

In the original self-taught learning framework [4|, the out-of-sample examples were used 
as a source of unsupervised data, and experiments showed its positive effects in a limited la- 
beled data scenario. However, many of the results by Raina et al. |4| (who used a shallow, 
sparse coding approach) suggest that the relative gain of self-taught learning vs ordinary su- 
pervised learning diminishes as the number of labeled examples increases. We note instead 
that, for deep architectures, our experiments show that such a positive effect is accomplished 
even in a scenario with a large number of labeled examples, i.e., here, the relative gain of 
self-taught learning is probably preserved in the asymptotic regime. 

Why would deep learners benefit more from the self-taught learning framework? The 
key idea is that the lower layers of the predictor compute a hierarchy of features that can be 
shared across tasks or across variants of the input distribution. A theoretical analysis of gen- 
eralization improvements due to sharing of intermediate features across tasks already points 
towards that explanation l24l . Intermediate features that can be used in different contexts 
can be estimated in a way that allows to share statistical strength. Features extracted through 
many levels are more likely to be more abstract (as the experiments in Goodfellow et al. Q 
suggest), increasing the likelihood that they would be useful for a larger array of tasks and in- 
put conditions. Therefore, we hypothesize that both depth and unsupervised pre-training play 
a part in explaining the advantages observed here, and future experiments could attempt at 
teasing apart these factors. And why would deep learners benefit from the self-taught learning 
scenarios even when the number of labeled examples is very large? We hypothesize that this 
is related to the hypotheses studied in Erhan et al. ||23l . Whereas in Erhan et al. |23| it was 
found that online learning on a huge dataset did not make the advantage of the deep learning 
bias vanish, a similar phenomenon may be happening here. We hypothesize that unsupervised 
pre-training of a deep hierarchy with self-taught learning initializes the model in the basin of 
attraction of supervised gradient descent that corresponds to better generalization. Further- 
more, such good basins of attraction are not discovered by pure supervised learning (with or 
without self-taught settings), and more labeled examples does not allow the model to go from 
the poorer basins of attraction discovered by the purely supervised shallow models to the kind 
of better basins associated with deep learning and self-taught learning. 

A Flash demo of the recognizer (where both the MLP and the SDA can be compared) can 
be executed on-line at |http : / /deep . host22 . com| 
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Appendix I: Detailed Numerical Results 



These tables correspond to Figures 2 and 3 and contain the raw error rates for each model and 
dataset considered. They also contain additional data such as test errors on P07 and standard 
errors. 

Table 1: Overall comparison of error rates (± std.err.) on 62 character classes (10 dig- 
its + 26 lower + 26 upper), except for last columns - digits only, between deep architec- 
ture with pre-training (SDA=Stacked Denoising Autoencoder) and ordinary shallow archi- 
tecture (MLP=Multi-Layer Perceptron). The models shown are all trained using perturbed 
data (NISTP or P07) and using a validation set to select hyper-parameters and other train- 
ing choices. {SDA,MLP}0 are trained on NIST, {SDA,MLP}1 are trained on NISTP, and 
{SDA,MLP}2 are trained on P07. The human error rate on digits is a lower bound because it 
does not count digits that were recognized as letters. For comparison, the results found in the 
literature on NIST digits classification using the same test set are included. 





NIST test 


NISTP test 


P07 test 


NIST test digits 


Humans 


18.2% ±.1% 


39.4%±.1% 


46.9%±.1% 


1.4% 


SDAO 


23.7% ±.14% 


65.2%±.34% 


97. 45 %±. 06% 


2.7% ±.14% 


SDA1 


17.1% ±.13% 


29.7%±.3% 


29.7%±.3% 


1.4% ±.1% 


SDA2 


18.7% ±.13% 


33.6%±.3% 


39.9%±.17% 


1.7% ±.1% 


MLPO 


24.2% ±.15% 


68.8%±.33% 


78.70%±.14% 


3.45% ±.15% 


MLP1 


23.0% ±.15% 


41.8%±.35% 


90.4%±.1% 


3.85% ±.16% 


MLP2 


24.3% ±.15% 


46.0%±.35% 


54.7%±.17% 


4.85% ±.18% 


H9| 








4.95% ±.18% 


|20| 








3.71% ±.16% 


[25] 








2.4% ±.13% 


[22| 








2.1% ±.12% 
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Table 2: Relative change in error rates due to the use of perturbed training data, either using 
NISTP, for the MLP1/SDA1 models, or using P07, for the MLP2/SDA2 models. A positive 
value indicates that training on the perturbed data helped for the given test set (the first 3 
columns on the 62-class tasks and the last one is on the clean 10-class digits). Clearly, the 
deep learning models did benefit more from perturbed training data, even when testing on 
clean data, whereas the MLP trained on perturbed data performed worse on the clean digits 
and about the same on the clean characters. 





NIST test 


NISTP test 


P07 test 


NIST test digits 


SDA0/SDA1-1 


38% 


84% 


228% 


93% 


SDA0/SDA2-1 


27% 


94% 


144% 


59% 


MLP0/MLP1-1 


5.2% 


65% 


-13% 


-10% 


MLP0/MLP2-1 


-0.4% 


49% 


44% 


-29% 



Table 3: Test error rates and relative change in error rates due to the use of a multi-task setting, 
i.e., training on each task in isolation vs training for all three tasks together, for MLPs vs 
SDAs. The SDA benefits much more from the multi-task setting. All experiments on only on 
the unperturbed NIST data, using validation error for model selection. Relative improvement 
is 1 - single-task error / multi-task error. 





single-task 


multi-task 


relative 




setting 


setting 


improvement 


MLP -digits 


3.77% 


3.99% 


5.6% 


MLP -lower 


17.4% 


16.8% 


-4.1% 


MLP -upper 


7.84% 


7.54% 


-3.6% 


SDA-digits 


2.6% 


3.56% 


27% 


SDA-lower 


12.3% 


14.4% 


15% 


SDA-upper 


5.93% 


6.78% 


13% 
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